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ABSTRACT 



Context. In the last years, the H 2 D + and D 2 H + molecules have gained great attention as probes of cold and depleted dense molecular 
cloud cores. These ions are at the basis of molecular deuterium fractionation, a common characteristic observed in star forming re- 
gions. H 2 D + is now routinely observed, but the search for its isotopologue D 2 H + is still difficult because of the high frequency of its 
ground para transition (692 GHz). 

Aims. We have observed molecular transitions of H 2 D + and D 2 H + in a cold prestellar core to characterize the roots of deuterium 
chemistry. 

Methods. Thanks to the sensitive multi-pixel CHAMP + receiver on the APEX telescope where the required excellent weather con- 
ditions are met, we not only successfully detect D 2 H + in the H-MM1 prestellar core located in the L1688 cloud, but also obtain 
information on the spatial extent of its emission. We also detect H 2 D + at 372 GHz in the same source. We analyse these detections 
using a non-LTE radiative transfer code and a state-of-the-art spin-dependent chemical model. 

Results. This observation is the first secure detection of D 2 H + in space. The emission is moreover extended over several pixels of 
the CHAMP + array, i.e. on a scale of at least 40", corresponding to ~4800AU. We derive column densities on the order of 10 12 - 
10 13 cur 2 for both molecules in the LTE approximation depending on the assumed temperature, and up to two orders of magnitude 
higher based on a non-LTE analysis. 

Conclusions. Our modeling suggests that the level of CO depletion must be extremely high (> 10, and even > 100 if the temperature of 
the core is around 10 K) at the core center, in contradiction with CO depletion levels directly measured in other cores. Observation of 
the H 2 D + spatial distribution and direct measurement of the CO depletion in H-MM1 will be essential to confirm if present chemical 
models investigating the basis of deuterium fractionation of molecules need to be revised. 

Key words, deuterium astrochemistry - prestellar core - depletion - submm observations 



1. Introduction 

Although deuterium is a trace element (D/H ~10 -5 ), the study 
of deuterium fractionation in molecules is a very vivid field of 
research. The main astrochemical interest derives from the fact 
that the understanding of fractionation processes may ultimately 
help to understand the formation of the related hydrogenated 
species. As an examp le, the fi r st de tection of a doubly-deuterated 
molecule, D2CO, by ITurnerl dl990h was interpreted as the hint 
that grain chemistry may be at work for synthesizing formalde- 
hyde. However, the impact of deuterium fractionation studies 
rapidly outpassed the mere astrochemical interest. Since deuter- 
ated molecules are formed preferentially at low temperatures and 
high densities, they became valuable probes for studying, e.g., 
the kinematics of cold cloud cores harbou ring the earliest stages 
of star formation (Ivan der Tak et al.ll2005l) . 

The basic process transferring deuterium from the HD 
molecule to other molecules is thro ugh ion-molecule reactions 
involving H+, CH+, and C 2 HJ (see IParise etal1l2009l and ref- 
erences therein). The route involving Ht is the dominant one 
at very low temperatures (T < 20 K), whereas the importance of 
CH3 takes over at slightly higher t emperatures (20 < T < 40 K) 
due to its higher endothermicity (ITurnerl 120011; iRouefF et alj 
I2007t IParise et al.1 120091) . H 2 D + is thus a very good probe of 
very cold objects, and has been extensively used in the last years 



* Based on observations with the APEX telescope. APEX is a col- 
laboration between the Max-Planck-Institut fur Radioastronomie, the 
European Southern Observatory, and the Onsala Space Observatory 



to probe gas in which CO is strongly deplete d, i.e., frozen out 
on dust grain surfaces (e.g. ICaselli et al.j 2003i), The deuteration 
level is also a decisive argument in the study of the evolutionary 
state of early phases of star formation. The deuterium fraction- 
ation is believed to de crease once the central pro tostar starts to 
heat its envelope (e. g., lEmprechtinger et al.l2009l). For instance , 
partly based on molecular deuteration. iBelloche et all {2006) 
concluded that the dense core Cha-MMS 1 in the Chamaeleon I 
molecular cloud is likely at the stage of the first hydrostatic core 
when H2 has not been dissociated yet. 

The study of several isotopologues of methanol toward the 
low-mass protostar IRAS 16293-2422 dParise et al.ll2002l |2004|) 
showed that grain chemistry models require a very high atomic 
D/H ratio accreting on the grains in order to explain the high 
CHD2OH and CD3OH abundances. At the time, gas-phase mod- 
els could not reproduce such a high atomic D/H ratio, but a 
decisive step forward was done with the inclusion o f D 2 H + 
and DJ in the reaction networks (iRoberts et al.l 12003b . Since 
then, D 2 H + was searched for through observation of its ground 
para state line. The frequency of this line w as initially pub- 
lished to be 691.660440 GHz (+ 19kHz (lcr\ iHirao & Amanol 
1 20031) and was later rev ised to 691.660483 GHz (+ 20kHz (lcr), 
Amano & Hiraol 12005b . These astronomical searches are ex- 
tremely difficult because of the high frequency of the molecular 
line, for which excellent instruments and weather conditions are 
required. Using t he Caltech Submill imeter Observatory 10.4m 
telescope (CSO). PVa"stel et all d2004l) detected a ~3.3cr line (in 
peak) toward the Ophiuchus dense core 16293E (4.4<x in inte- 
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grated intensity). This detection is to date the only published 
one, and its low signal-to-noise ratio as well as the velocity shift 
between H2D + and D2H + (0.23 kms" 1 after remeasure ment of o- 
H2D + and p-D2H + frequencies, Amano & Hiraoll2005l) certainly 
require a confirmation of the detection by improving the ob- 
servation sensitivity and frequency calibrat ion precisio n (whic h 
was ~0.1kms _1 in the observations of Vast el et a D (12004)). 
We present here the first secure detection of D2H + toward the 
prestellar core H-MM1 in L1688, the main molecular cloud in 
the Ophiuchus star forming region. 

2. Observations 

2.1. The source 




RA(2000) 



Fig. 1. The L1688 region, as mapped by SCUBA (850 pm), from 
the SCUBA legacy catalogues (see text). The contours are in 
steps of 3<x (=0. 18 Jy/beam). The beamsize is 14" (HPBW). 



The L1688 cloud (the main cloud in the Ophiuchus 
dark cloud) has b een mapped both in molecular tracers (e.g. 
lLoren et aL 19901) an d in continuum emission (e.g. at 1.3 mm 
bv iMotte etal] 11998b . These early mapping studies revealed 
different areas of high column density, named Oph-A to 
Oph-G. Using the Submillimeter Common User Bolometer 
Array (S CUBA) at the 1 5 m Jam es-Clerck-Maxwell Telescope 
(JCMT), iJohnstone et alj d2004l) covered a larger region, lead- 
ing to the detection of two new dense cores, that they named 
Oph-H and Oph-I, following the same nomenclature. Oph-H 
contains a single continuum peak, H-MM1. The distance of the 
Ophiuc hus region has been recently accurately determined to 
120pc dLoinard et alJl2008h . 

We present in Fig. Q] a part of the SCUBA map 
of the Ophiuchus region, as available from the SCUBA 
Legacy catalogued ] (data: scuba_F_353dl_16d8_850um.emi.fits, 
iDi Francesco et al" 120081) . showing the relative position of 
all dense cores. In this paper, we present observations to- 
ward the H-MM1 core, at coordinates ff 2 ooo=16' ! 27'"58.3 s and 
S 2 ooo=-24°33'42.0". 

H-MM1 corresponds to the Bolo2 5 source of the 1 . 1 mm 
BOLOCAM (CSO) map observed bv lYoung et al] d2006l) . The 



FWHM beam size of these observations is 31". Assuming a 
temperature of 10 K and a dust opacitjQ k v = 0.0114 cm 2 g _1 
at 1.12mm these authors infer a total mass of O.24M in 40", 
and measure a source FWHM size of 47"x56". They measure 
a peak visual extinction, Ay, of 23 mag, and a mean density of 
5xl0 5 crrT 3 . 



iKhanzadvan et all d2004l) and lStanke et al.1 d__006) carried out 
an unbiased search for protostellar signatures in the p-Oph cloud. 
The closest source they identify in their continuum maps is 
MMS38 (at 9" from our position), which they classify as star- 
less. No optical/NIR counterpart (2MASS) or sign of outflow 
are found in this source. 

The region was also mapped with the Spitzer Space tele- 
scope, in the frame of the Cores to Disks legacy project (c2d). 
We downloaded the corresponding data from the Spitzer web- 
siteS No source is detected in any of the Spitzer bands, espe- 
cially _a^70jwrvwhichjs a good tracer of embedded protostars 
(see iDunham et al1l2008h . H-MM1 thus appears to be a dense 
prestellar core. 

2.2. APEX line observations 

Using the APEX2a single-pixel and CHAMP + array receivers at 
the Atacama Pathfinder Experiment telescope (APEX), we ob- 
served H 2 D + (372 GHz) and D 2 H + (692 GHz) toward H-MM1 
at the position given in Sect. 2.1. We used in each observation the 
position-switch mode, with the off-position at (0", 300") equa- 
torial offsets. Calibration was done by measuring regularly the 
sky brightness and cold and hot loads in the cabin. 

2.2.1. H,D + 



Using the APEX2a receiver (Risac her et al.l [2006), we ob- 
served the ground-state o-H2D + line to ward H-MM1. The fre - 
quency of this line is 372.421385 GHz dAmano & Hiraoll2005l) . 
Observations took place on June 29th, 2007, under very good 
weather conditions (precipitable water vapor PWV=0.24mm). 
The system temperature was 250 K. The backend, a Fast Fourier 
Transform Spectrometer (FFTS), was used with 8192 channels, 
yielding a 0.0983 km s _1 channel separation. The effective fre- 
quency res olution of the FFTS in this configuration is 159 kHz 
dKlein et alJl2006l and B. Klein, priv. comm.), i.e. 0.13km/s at 
372 GHz. 

The beam efficiency of APEX at 372 GHz was measured to 
be 0.73 for compact sources (as measured on Mar s and Jupiter, 
which were respectively 8" and 32-35" at that time. lGusten et al.l 
2006). As the D2H + emission seems extended on a scale compa- 
rable to the footprint of CHAMP+ (-40", see Figure 0J, H 2 D + 
is likely to also be extended on that scale. We therefore adopt 
the T„,t, scale for our analysis, using the beam efficiency of 0.73. 
The forward efficiency of the telescope is 0.95 at this frequency. 

We reached a rms noise level of 72 mK (T* scale) at the ini- 
tial resolution (0.13kms _1 ). 

2.2.2. D 2 H + 

Using the CHAMP + array receiver on the APEX telescope, we 
targeted the prestellar core H-MM1 . The observations took place 



1 http://www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/ 
scubalegacy/ 



2 This dust opacity assumes a gas-to-dust mass ratio of 100, and is 
expressed per dust+gas mass unit 

3 http://data.spitzer.caltech.edu/popular/c2d/20071 101_enhanced_vl/ 
oph/ 
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RA offset (arcsec) 



Fig. 2. 850 /im continuum emission of the H-MM1 core. The 
map is centered on the position of H-MM1 (see text). The con- 
tours are in steps of 3<x(=0.18 Jy/beam). The blue numbers show 
the different pixels of the CHAMP + low-frequency array. The 
blue circles show the D2H + beams. In red, the H2D + beam. 

in June and August 2009, under good to very good weather con- 
ditions (most of the time with PWV < 0.45mm). The typical sys- 
tem temperature in the central pixel of the array varied between 
800 and 1400 K. Figure presents the overlay of the CHAMP+ 
footprint on the 850 /im continuum emission. Each CHAMP + 
pixel was tuned to the CO(6-5) frequency (691.4730763 GHz), 
which is separated by 187.4MHz from the D 2 H + frequency. 
We ad opt here the most recent measurement of Ama no & Hiraol 
(120051) of the D 2 H + line frequency: 691.660483 GHz. The beam 
size is 9". The CHAMP + AFFTS was used with a channel sep- 
aration of 183 kHz (0.079 km/s). The effective resolution is 
212kHz, i.e. 0.092 km/s. 

The beam efficiency at 692 GHz was measured to be 0.52 
on Jupiter (angular size 47.3" at that time). Because the D 2 H + 
emission is observed to be extended (see below) on a typical 
size of 40", we adopt the beam efficiency of 0.52. The forward 
efficiency is 0.95 at this frequency. 

We reached an rms noise level of 34 mK (T* scale) at the 
initial resolution (0.09 km s -1 ). 




Velocity (km/s) 



Fig. 3. H 2 D + (upper panel) and D 2 H + (lower panel) toward the 
central position of H-MM1. 




2.3. Continuum data 



We 



retrieved the 850/im map of the Ophiuchus re- 
gion from the SCUBA Legacy catalogues (data: 
scuba_F_353dl_16d8_850um.emi.fits iDi F rancesc o et alj 

120081) . The region co vering the L1688 cloud is shown in Fig. [T] 
and was discussed bv lj0rgensen et alJ ([2008). 

The beam of SCUBA at this frequency is 14". The local rms 
in the H-MM1 region of the map is found to be 60 mJy/beam. 

3. Results 

The main result of this work is the first unambiguous detection of 
p-D 2 H + toward a prestellar core. Figure[3]presents the H 2 D + and 
D 2 H + emission at the (0", 0") position. The agreement of the 



Fig. 4. D 2 H + observations. Each panel shows a pixel of 
CHAMP + (as displayed on Fig. [2]). The vertical line shows 
the vis, of the source as measured with the o-H 2 D + line 
(4.23 km s~'). The y-axis is in Kelvin, on the T* scale. 



LSR velocities, v/ jr , for both lines is remarkable (see Table Q]) 
and further strengthens the assignment of the lines. Figure [4] 
presents the D 2 H + observation on each of the CHAMP + pix- 
els. The noise level is different from pixel to pixel, due to the 
different performances of the off-axis pixels. The line is clearly 
detected in at least two adjacent pixels (pixels 2 and 4). There is 
also a tentative detection in pixel 3, with a lower intensity. The 
D 2 H + emission is thus extended on the scale of the CHAMP + 
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Table 1. Observational results. 



Species 


Position 


Pixel 


\ T* dv 




Av 


a, peak 


[A v](/ econv 


5* 


I T* dv / 5 




EQ offsets 




(mK km s ) 


(km s~') 


(km s~') 


(K) 


(km s~') 


(mjy/bea.ni) 


(K km s -1 / [Jy/beam]) 


o-H 2 D + 


(0",0") 


1 


180±20 


4.23+0.04 


0.51+0.07 


0.33 


0.49 






p-D 2 H + 


(0",0") 


1 


47±8 


4.30+0.03 


0.35+0.06 


0.13 


0.34 


628+60 


0.076+0.015 




(9.5", -17") 


2 


49±8 


4.26+0.04 


0.42+0.09 


0.11 


0.41 


468+60 


0.106+0.022 




(19.5", -1") 


3 


27+10 


4.16+0.05 


0.31+0.10 


0.08 


0.30 


404+60 


0.068+0.027 




(10", 16") 


4 


44+11 


4.23+0.05 


0.41+0.12 


0.10 


0.40 


616+60 


0.072+0.019 




(-9.5", 17.0") 


5 


< 30 










343+60 


< 0.087 




(-19.5", 1") 


6 


<20 










211+60 


< 0.095 




(-10", -16.3") 


7 


<20 










179+60 


< 0.112 



* S is the peak flux density in the 850/mi SCUBA map at the center of each CHAMP + pixel. It is in units of mJy/14"-JCMT beam. The noise in 
the map is 60 mjy/beam. The error bars are all lcr, and the upper limits 3cr. 



footprint (~40"), corresponding to a physical size of ~4800 AU 
at the distance of 120 pc. The observational results are listed in 
Tabled] 

We note that the detection of p-D 2 H + in H-MM1 does not 
imply that it is a peculiar core. No ne of the searches for p-D2H + 
published in Caselli et al. (2008) towards five prestellar cores 
have reached upper limits which are deep enough to discard the 
detection of the line with the same peak temperature in those 
cores. 



3.1. Linewidths 

The measured linewidths are listed in Table Q] Because the lines 
are not very broad compared to the frequency resolution of the 
FFTS, we also computed deconvolved linewidths to correct for 
instrumental broadening (column 7). 

The expected thermal linewidth for a molecule of mass m is: 



AV; 



therm.FWHM 



2 x V2 ln2 x -\ / — 

m 



For H 2 D + , Av therm = 0.107 x Vrkm/s. For D 2 H + , 
Avtherm = 0.0955 x Vf km/s, with T in Kelvin. 

The small linewidth of D 2 H + on the (0", 0") position con- 
strains the kinetic temperature to T < 1 3 K. For this upper limit 
value of the temperature, the rather large H 2 D + linewidth points 
to non-thermal velocity broadening of 0.3 km/s. The different 
non-thermal velocity contribution to both lines might be a sign 
that the two molecules do not trace exactly the same gas. The 
larger non-thermal velocity broadening of o-H 2 D + may imply 
that o-H 2 D + is present in a more diffuse and turbulent gas than 
p-D 2 H + , as also supported by chemical models, which predict 
that the p-D 2 H + / o-H 2 D + ratio raises with density (see Figure|9]). 
In this case, the p-D 2 H + / o-H 2 D + ratio that we will derive in the 
densest gas will likely be a lower limit. 

3.2. Correlation with continuum emission 

The continuum emission at the position of the individual pixels 
of CHAMP+ is listed in Table Q] The positions where p-D 2 H + is 
detected correspond to the positions where the continuum emis- 
sion is the strongest. The lower tentative detection in pixel 3 is 
also consistent with the continuum being weaker at this posi- 
tion. The correlation between the p-D 2 H + and the continuum 
emission is quite good. This implies that the p-D 2 H + emission 
comes from the densest regions of the core, in agreement with 
the narrow width of the line. 



3.3. Column densities 
3.3.1. LTE assumption 

The critical densities for the o-H 2 D + and p-D 2 H + transition s 
are, respectively, 1.3x10 s and 5.6xl0 5 ctrT 3 dHugo et al.ll2 009). 
Although the density of the dense core may be of the same order 
or only marginally higher than these critical densities, we start 
with an LTE analysis as a first approximation. We will consider 
a non-LTE approach in the following section. 

We consider the ortho and para species of each molecule as 
different species, that do not interconvert via collisions with H 2 . 
We thus compute separately the partition function of ortho and 
para species, from the transitions given in the Cologne Dat abase 
of Molecular Spectroscopy (CDMS, iMiiller et alJ l2005). The 
partition functions are computed as 



Q(T) = J] 



gig r e 



(-E/kT) 



levels 

where gj is the spin degeneracy (1 and 3 for para- and ortho- 
H 2 D + respectively, and 1 and 2 for para- and ortho- D 2 H + ) and 
g r the rotational degeneracy (g r — 2J + 1). Note that the ener- 
gies of the levels were all expressed relative to the ground level 
Ooo of each molecule. These partition functions are presented in 
Table |2] 

Table 2. Partition functions of o-H 2 D + and p-D 2 H + at different 
temperatures. 



Temperature 


Q(o-H 2 D + ) 


Q(p-D 2 H + ) 


(K) 






7 


4.2xl0~ 5 


2.3xl0~ 3 


8 


2.0xl0~ 4 


5.7xl0" 3 


9 


6.9xl0~ 4 


1.2x10^ 


10 


1.8xl0~ 3 


2.0xl0~ 2 


11 


4.2xl0~ 3 


3.3xl0 2 


12 


8.2xl0~ 3 


4.8xl0~ 2 


13 


1.5X10 2 


6.8xl0~ 2 



The column density in the upper level of the transition can 
be directly derived from the measure of the opacity. 



N u = 



1 



1 



hB u i e hv ' kT - - 1 



» peak 



Av 
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Table 3. Line opacities and column densities on the (0", 0") position in the LTE hypothesis. 



T 

1 ex 


t(o-H 2 D + ) 


jV(o-H 2 D + ) 


r(p-D 2 H + ) 


N(p-D 2 H + ) 


yV(p-D 2 H + )/^V(o-H 2 D + ) 


(K) 




(cm" 2 ) 




(cm" 2 ) 




7 


0.35 


9.3xl0 12 


1.7 


4.1xl0 13 


4.4 


8 


0.23 


6.5xl0 12 


0.59 


1.5xlO B 


2.2 


9 


0.17 


5.1xl0 12 


0.33 


8.2xl0 12 


1.6 


10 


0.13 


4.2xl0 12 


0.21 


5.4xl0 12 


1.3 


11 


0.11 


3.6xl0 12 


0.15 


3.9xl0 12 


1.1 


12 


0.088 


3.2xl0 12 


0.11 


3.0xl0 12 


0.96 


13 


0.076 


2.9xl0 12 


0.089 


2.5xl0 12 


0.86 



The opacity of the transition, r, is given by: 



' mb 



[J(T ex )-J(T bg )][l-e- T ] 



where J(T) 



hv/k 



exp(hv/kT)-l 

Assuming LTE (i.e. T ex = 7\i n ), 
column density of the species: 



(1) 



we can then derive the total 



Ntn, = 



EJkT„ 



' peak 



Av 



guA u i e hv ' kT - - 1 

We used the A u / of iRamanlal & Tennyson! (2004). For both 
transitions, they are 1 1 % different from the value quoted in the 
CD MS database . We u sed the accurate frequencies derived by 
lAmano&Hiraol (120051) . 

We note that equation (1) implies that r m b<J(r ex ). This 
condition for the 692 GHz line of D2H + translates into 
J(7ex) ^ 0.25 K, i.e. T ex > 7 K. This condition as well as the pre- 
vious linewidth argument allow to restrict the kinetic tempera- 
ture to the range 7-13 K. 

The results are presented in Table [3] Depending on the as- 
sumed temperature, the P-D2FF/ o-H2D + ratio varies in the 
range 0.86 - 4.4. 

3.3.2. Out-of-equilibrium radiative transfer 

Thanks to the recent comput ation of all elastic and non-elastic 
collision rate coefficients bv iHugo et alJ ((2009), we can derive 
the column density using an out-of-equilibrium radiative trans- 
fer model. We used the n on-LTE radiative transfer code RAD EX 
(Ivan der Tak et alJ [2007), in the isothermal sphere geometry, 
with the escape probability as a function of the opacity as fol- 
lows: 



1.5. 



-II 



- + (- + -)e- T ] 

T L 7- 7" z 



Because the rates for collisions with 0-H2 and P-H2 are quite 
different, the nature of the collider matters for the excitation. Our 
chemical models (presented in Sect. [U predict that the steady- 
state ortho/para-H2 is within a factor of two around 10~ 4 for 
all physical conditions probed (Fig. [8}. We checked that the ex- 
citation is unchanged (within numerical precision) when con- 
sidering only para-H2 or para-H2 with traces of ortho-H2 (o/p- 
H2<2xl0~ 4 ). We thus assumed in a first approximation colli- 
sions with para-H2 only in the derivat i on of the column density 
of o-H 2 D + and p-D 2 H + . Pag ani et all J2009h however showed 
that the 0/P-H2 ratio is unlikely in steady-state in prestellar cores 



p-D 2 H + □ n = 10 4 cm" 
o-H 2 D + a n = 1 5 cm" 
o n = 1 5 cm" 




9 10 11 

Kinetic temperature (K) 



i / 



- b) 













7 8 9 10 11 12 13 

Kinetic temperature (K) 

Fig. 5. a) Column densities at the (0,0) position, computed from 
the non-equilibrium radiative transfer assuming the beam fill- 
ing factor for both molecules is 1, as a function of the physical 
conditions (kinetic temperature, density). The thick black lines 
show the column densities computed in the LTE assumption, b) 
Excitation temperature in the same conditions and color conven- 
tions as Figure a. The black line represents T ex = 7\j n . 



(see also Section l4~2l and Fig. [7]i- F° r completeness, we checked 
the computed column densities in the case where the 0/P-H2 ra- 
tio is higher than the steady-state value. We find that the derived 
column densities start to deviate sensibly from the ones derived 
with only P-H2 for values of the 0/P-H2 ratio of the order of unity. 

The results with excitation with para-H2 only are presented 
in Fig. [5^. The results with an 0/P-H2 ratio of 10 -2 are found to 
differ from Fig. [5^ by less than 1%. The column densities derived 
within the LTE approximation are also shown on the figure. Even 
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at a density of 10 6 cm 4 , the LTE assumption is not accurate, 
especially for p-D2H + which has the highest critical density. One 
can check indeed that the excitation temperature does not reach 
the kinetic temperature in any of the studied cases (Fig. [5J3). 

Assuming 0/P-H2 = 3 (the highest possible ratio at thermody- 
namical equilibrium) results in column densities about 10-20% 
higher than in the case of excitation with only para-H2, the ef- 
fect being more important at low densities, as expected (at high 
densities, the system is close to LTE, and does not depend on the 
nature of the collider). The abundance ratio of the two molecules 
is even less sensitive to the 0/P-H2 ratio, as shown in Fig. [6] This 
figure presents the p-D2H + / o-H2D + ratio resulting from our ra- 
diative transfer analysis, for different densities and in the two 
cases 0/P-H2 = and 3. In the following, we aim at investigating 
the implications of this the p-D2H + / o-H2D + ratio with the help 
of a state-of-the art chemical model. 
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Here hq and crc are the density and cross section of grain par- 
ticles, Vh, «h, Vd, and «d are the thermal velocity and number 
density of H and D atoms. In this limit, the rather controversial 
detailed formation mechanism of H2 and the uncertain evapora- 
tion and diffusion energy barriers do not enter the equations. The 
ortho and para forms of H2 and D2 are produced according to the 
statistical weights (3 : 1 and 2:1). 

The influence of the different pa rameters on the chem - 
istry has been studied in great de t ail by Walmsl ev et al] d2004). 
iFlower et al.l(l2004l) . lPagani et al.l(l2009l) and lSipila etafl hold) . 
Here, we adopt a cosmic ray ionization of 3xl0~ 17 s _I , a dust-to- 
gas mass ratio of 1.3X10 -2 , and a grain radius of 0.1 fim. These 
parameters are not varied, because of the lack of observational 
constraints at hand. 

We simulate different levels of CO depletion by using dif- 
ferent CO initial abundances. Note that we run the models with 
a constant elemental C and O abundance (i.e. we do not allow 
depletion on the grains during the evolution). In all our models, 
we assume that depletion affects CO and N2 abundances at the 
same level. 



(k) 



Fig. 6. p-D2H + / o-H2D + abundance ratio resulting from the radia- 
tive transfer modelling of the observations, for different densities 
and 0/P-H2 ratio. 



4. Astrochemical Modelling 

4.1. Our model 

We solve the chemistry in the gas phase using a code devel- 
oped by our group, and based on the DLSODES solver, which 
is part of the ODEPACK package (w ww.netlib.org) . The code 



was benchmarked against the results of iPagani et al. 



(200 91 and 



Sipil a et aTl d2010l) . We used the same approach as Sipil a et al.l 
d2010l) in building the reaction file for the complete depletion 
case. We take into account all different spin states of the different 
molecules and ions (i.e. we consider separately ortho and para 
species of H2, D2, H2D + , D2H + , H^ and ortho, para and meta 
sp ecies of ) . We use the reaction rates between these species 
of lHugo et all (f2009). In addition we a dded a simple chemistry 



of CO and N2, as in Pagani et al. ( 20091). W e use the recombina 



Paga ni et all (12009) 



tion rates for all ions of 

The pseudo-time-dependent chemical model contains only 
gas phase reactions, except that formation of H2, HD, and D2 
on the grain surface is als o included. We use the "large grain 
limit" (Lipshtat et al. 2004), with the further assumption that all 
the influx of H and D atoms is converted into H2, HD, and D2 
molecules, and the diffusion rates of H and D atoms on the grain 
surface are approximately the same. With these assumptions, the 
formation rates of these three molecules can be written explicitly 



4.2. Timescales 

The free-fall timescale T/f of an isothermal sphere of 
n = 10 6 cirT 3 is 3xl0 4 yrs. Magnetic fields can increase the dy- 
namical timescale, then controlled by ambipolar diffusion in the 
early evolution of subcritical cores. Timescales of core evolution 
can be inferred empirically by comparing the number of prestel- 
lar and protostellar cores in observed regions (for a review, see 
IWard-Thompson et al.l2007l and references therein). These em- 
pirical timescales lie roughly in the range 2-5 x tjj. 

It is instructive to compare the chemical timescales to this 
dynamical timescale. Figure [7] shows the time evolution of the 
0/P-H2 ratio and p-D2H + /o-H2D + ratio for constant physical 
conditions. The temperature was fixed to 12K, and three dif- 
ferent densit i es wer e considered. As already pointed out by e.g. 
Pagani et al. (2009), the longest timescale in the chemistry is 
the conversion between the two forms of ortho and para H2. 
The other ions adjust to the o/p ratio of H2 in typically less 
than 10 2 yrs. The ortho-para interconversion of H2 takes around 
10 4 yrs (typical timescale for decrease of 0/P-H2 ratio by a factor 
1/e, see Fig.[7j). 

The free-fall timescale varies as l/V" an d the chemical 
timescale varies even slower with the density. The timescale 
for 0/P-H2 conversion is 6xl0 3 yrs at density n — 10 6 cm~ 3 , and 
9xl0 3 and 1.7xl0 4 yrs for densities of 10 5 and 10 4 cm -3 , respec- 
tively. At high densities, the chemical timescales are somewhat 
shorter than the free-fall timescales (thus roughly one order of 
magnitude smaller than the empirical core timescales), and the 
difference even increases at lower densities. This comparison be- 
tween timescales justifies at first order to consider the chemistry 
in dependently of the dyn amical evolution, as already discussed 
bv lWalmslev et all (120041) . 

Because the evolutionary stage of the H-MM1 core is not 
constrained, in the following we will focus on the study of the 
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Fig. 7. Time evolution of the chemistry in the case of T = 12 K. 
The different colors stand for different densities: n=10 4 (blue), 
10 5 (green), 10 6 cm ~ 3 (red). The thick lines represent the o/p 
H2 ratio. The thin lines represent the p-D2H + /o-H2D + ratio. The 
different line styles stand for different CO abundances: complete 
depletion (full), x co = 10~ 6 (dash), x co = 10~ 5 (dash-dot). 
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Fig. 8. Steady-state ortho-to-para ratio of H2. 



chemical composition when the steady-state is obtained. This 
choice avoids including the age of the core as an additional pa- 
rameter but will not invalidate the generality of our conclusions. 
As shown in Fig. [7j taking the chemical composition at long 
times overestimates the p-D2H + /o-H2D + ratio in case the core is 
not yet in chemical steady-state. This note will be important for 
the extrapolation of our conclusions in the next sections. 

4.3. Physical conditions 

We could in principle derive the density profile of the core 
from the continuum emission at 850//m. However, this would 
require the knowledge of the te mperature profile, w hich is not 
known. The Gould Belt Survey (And re et alj|2010b made with 
the Herschel Space Observatory should allow to derive accurate 
temperature profiles. Until then, as we cannot constrain the den- 
sity profile without ad-hoc assumptions on the temperature, we 
restrict ourselves here in running chemical models at different 




7 8 9 10 11 12 13 

Kinetic temperature (K) 



Fig. 9. Chemical model predictions of the p-D2H + /o-H2D + ratio 
(same convention as previous figures), at steady-state. The thick 
decreasing lines show the ratio derived from the non-LTE anal- 
ysis of the observations, in the assumption of different densities 
and different 0/P-H2 ratios (same curves as presented in Fig. [6] 
but with the same linestyle for both o-p H2 ratios). 



densities, temperatures, and levels of CO depletion. Our goal is 
to investigate under which average conditions the model can re- 
produce the observations. A full physical model of the source 
will be presented in a forthcoming detailed study, after collect- 
ing new observational constraints. 

Figure|9]shows the prediction of the model for the p-D2H + /o- 
H2D + ratio. The thin lines are the model predictions, as a func- 
tion of temperature, and for different densities and CO depletion 
factors. At fixed density, the p-D2H + /o-H2D + ratio is increas- 
ing with the CO depletion level, because HD becomes more and 
more competitive over CO for reaction with and isotopo- 
logues. The ratio increases with tempe rature in this t emper ature 
range, as already explained in detail bv lFlower et al . (2004- The 
three thick lines are the ratios as calculated with RADEX from 
the observed line intensities (as a function of temperature and for 
different densities and 0/P-H2 ratios, as already shown in Fig.|6]l. 
We discuss in the following section the comparison between the 
chemical model and the observations. 



5. Discussion 

Figure [9] allows to investigate the effect of the temperature, den- 
sity and CO depletion level on the P-D2FF / o-H2D + ratio, and 
to infer under which average conditions our model could repro- 
duce the observed value. It should be noted that the observa- 
tional ratio also depends on the assumed density, because of the 
non-LTE effects described in Section [3.3.2l Agreement between 
observations and model can only be reached for densities strictly 
higher than 10 5 cm 3 . This is consistent with the re latively high 
measu red average density of the core (5xl0 5 cm~ 3 . iYoung et alJ 
2006). Even at the high density of 10 6 cm 3 , the model can only 
reproduce the large observed p-D2H + /o-H2D + ratio if T >10K 
and if the CO depletion level is substantial. The model points 
to the fact that the CO abundance should be < 10~ 5 . This corre- 
sponds to a CO depletion level of more than 10 (for T — 13 K). 
The required depletion is even more severe if the density is a bit 
lower, or if the temperature is closer to 10 K. Complete depletion 
(i.e. »100) is required at 10K. 
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iBacmann et al.l (2002) have studied the level of CO deple- 
tion toward a sample of seven prestellar cores. These cores are 
located at distances between 120 and 200 pc, and three of them 
are in the Ophiuchus molecular cloud. The cores have mean den- 
sities in the range (0.9-6) x 10 5 cm 4 , similar to H-MM1. Based 
on IRAM observations of dust continuum and CO isotopologues 
(at angular resolution of 22"), the authors find depletion factors 
at the continuum peak of the cores in the range 4.5-15.5, but 
this factor drops to 2-5 at offset positions where the continuum 
is 70% of the peak value. Although the study of Bacmann et 
al. (2002) was done at a 22" resolution (to be compared to our 
9-14" beams), preventing us from doing an accurate compari- 
son, the average depletion level we deduce in H-MM1 seems to 
be in the high range. Moreover the D2H + emission is extended 
over 40". Should the H2D + be also extended on that scale, this 
would imply that the model requires a high CO depletion (>> 
10 ) over an extended re gion, in contradiction with the results 
of Bacma nn et alj d2002l) . It is clear that the observation of the 
H2D + spatial distribution and the direct measurement of the level 
of CO depletion are required to confirm if the astrochemical 
model is indeed unable to describe satisfactorily the roots of deu- 
terium chemistry. 

Following the discussion on the linewidths in Section lXTl the 
level of the thick observational curves in Fig. [9] is in fact likely 
underestimated. Similarly, following the discussion in Section 
14.21 the model predictions are overestimated if the core has not 
yet reached chemical equilibrium. These remarks show that the 
difficulties of the model might be even more severe than shown 
on Fig. [9] 

The model also seems to exclude kinetic temperatures be- 
low 10K, but several v e ry cold prestellar cores are known (6 - 
7 K, ICransi et alJ l2007t iPagani et alJ l2007t lHariu et al.l [2008). 
The temperature in H-MM1 still needs to be determined, but this 
point could be another difficulty of the model. 

Further tests of our understanding of the root of deuterium 
fractionation will depend on the observation of H2D + at the off- 
set positions where D2FF was detected, and on the independent 
measure of the kinetic temperature and CO depletion level, for 
which we have an on-going APEX project. 

6. Conclusions 

We presented the first secure D2H + detection in the ISM, to- 
ward the Oph H-MM1 core. The centroid velocity of the line is 
consistent with the one of the o-H2D + line. The emission is ex- 
tended over 40" at least. Using an out-of-equilibrium radiative 
transfer model, we have inferred the column densities of both 
o-H2D + and p-D2H + . We analyzed the chemistry with a state- 
of-the-art astrochemical model including all spin states of , 
H2D + , D2I-F, D3 . The model can reproduce the data only if the 
average density is high (n > few 10 5 cm ~ 3 ), the kinetic temper- 
ature is between 10 and 13 K, and the depletion level of CO is 
high (> 10 at 13 K, > 100 at 10K). Further observational studies 
may show whether these conditions indeed apply to H-MM1, 
and if the CO depletion level stays high on the large scales at 
which p-D2H + emits. 

Should new observational studies confirm the failure of the 
model, new ways will have to be investigated. We list here sev- 
eral possibilities: include a more complete description of the 
chemistry (treating fully e.g. the CO chemistry and gas-grain in- 
teraction), recompute or measure the reaction and collision rates 
of all (spin separated) species, include in the chemical model a 
correct treatment of the excitation of the different molecules (at 



the moment, all adopted reactive and non-reactive collision rates 
describe the molecules in their rotational ground state). 
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